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We investigate the boson star with the self-interacting scalar field as a model 
of galactic halos. The model has slightly increasing rotation curves and allows 
wider ranges of the mass(m) and coupling(A) of the halo dark matter particle 
than the non-interacting model previously suggested(ref.|0|). Two quantities are 
related by \^{mp/m)'^ ~ 10^°. 







It is well known that the flatness of the galactic rotation curves indicates the 
presence of dark matter around galactic halos. However, the properties of the 
dark matter are still mysterious. For example, why the dark matter in halos does 
not fall towards the center of galaxy and form black holes? The answer to the 
above question may be a good criterion for the good halo model. 

There are thermal distribution model [|l| where density profile p ~ r^^ , and 
spherical infall model where p ~ r^^-^^. 

Recently Sin[^ Q suggested a new model of the halos composed of pseudo 
Nambu-Goldstone boson (PNGB). 

According to the model, the condensation of ultra light PNGB whose Comp- 
ton wavelength \comp = is about Rhaio is responsible for the halo formation. 
Cosmological role of the ultra light PNGB was studied in the late time phase tran- 
sition model to reconcile the smoothness in the background radiation with the 
large scale structure. 

Before Sin's work, an astronomical object which consists of the PNGB dark 
matter was suggested by some authors]^. In their model the force against gravita- 
tional collapse comes from the momentum uncertainty of the quantum mechanical 
uncertainty principle. 

Since the typical length scale R in this model is Compton wavelength \comp ~ 

2 

^ of the particle, the typical mass scale of the object is M ~ ^ ~ 

Similarly, in Sin's model galactic halos are the objects of the self-gravitating 
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bose liquid whose collapse are prevented by the uncertainty principle. 

The typical halo has radius Rhaio ~ lOOkpc ~ lO^^cm and mass Mhaio ~ 
~ ^O'^^g, so one find the mass m of the PNGB whose de Broglie wave 
length ~ Rhaio is about 10~'^^eV. 

Note that the de Broglie length ~ f x Xcomp is more adequate to our purpose. 

The self-gravitating condensed states are described by the following non-linear 
Schroedinger equation: 

ihdtiJ = -— V^^ + GmMo r dr'^ f dr"A-Kr"'^\i)\^i)ir), (1) 
2m Jo r'^ Jo 

which was known as the Newtonian limit of the boson star fields equation]?]]. 

The normalization constant Mq is chosen to give the total mass of halo M = 
Mo / drAirr^l^pl'^ as in ref.|||. 

The rotation velocity of the stellar object rounding halo at radius r is given 

by 

^ (2, 

where M(r) is mass within r. 

Integrating eq.(|l]) numerically and using eq.(||) Sin found slightly increasing 
rotation curves and density profile p ~ r~^'^. 

What happens if there are repulsive self-interactions between the dark mat- 
ter particles? To answer this and stability question it is desirable to study the 
relativistic fields equations than the Schroedinger equation. 
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The cold gravitational equilibrium configurations of massive scalar field were 
found by solving the Klein-Gordon equations with gravity decades ago[|]. We find 
that these configurations, called boson star 0, are adequate to the relativistic 
extension of Sin's model. 

Consider a self-interacting complex scalar field and the gravity whose action 
is given by 



Since halos seem to be spherical, we choose Schwarzschild metric 

ds^ = -B{r)df + A{r)dr^ + r^dVL (4) 
and assume spherically symmetric field solutions 

0(r,t) = (47rG)-^(T(r)e-*"*. (5) 

From the action, dimensionless time independent Einstein and scalar wave equa- 
tions appear as in ref . [|10(] : 

A' 
B' 









(6) 
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B' 
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\a' + A[i^-l)a-Aa']=0, 


(8) 



where x = mr, n = — , A = \1 — 2^^^^^] ^ and A — -. — y • 

One may take M{x) for dimensionless mass of the boson star for large x. 
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Numerical solutions of the above equations are studied by many authors |TT], 



T2| , p!3| . The required boundary conditions are M(0) = 0, cr'(O) = and B{oo) 



1 and free parameters are cr(0) and Q. 



For the case A = |jT2[ it was found that there is a maximum mass Mj^ax = 
0.633— for the zero node solution. 

m 

We will focus on the non-zero node solutions, because the rotation curve of the 
zero node solution falls too fast to explain the flatness of the rotation curves of 
many galaxies(see ref. |^ for more arguments). This raises the stability problem 
of higher node solutions, which will be discussed later. 

Maximum masses for higher node solutions are proportional to node number 
n and about the same order as for the zero node case for small n. 

This with Mhaio gives us m ~ 10~'^'^eV. 

Another constraint comes from the maximum stable center density against 
small radial perturbation ||13| pc = 2.1 x lO^^m^g/cm^ > IQ^^^g/crn? , which is 
equivalent to m ~ 10~^^el^ for the zero node solutions. So for the zero node 
solutions 10-28eF ~ m ~ l{)-'^'^eV. 

For the case A 7^ 0, new scale appears because of the repulsive force preventing 
halo from gravitational collapse. In this case the typical length scale is i? ~ A^/m, 
thus the typical mass scale is ^ ~ A2 m^/m, which is also of order maximum mass 
like A = case. 

Numerical study||lO| shows M„iax = 0.22A3^ for zero node solutions. From 
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the fact that M^ax > Mhaio we find 

A^(!!^)2 > 1050. (9) 
m 

This is a relation between the mass and coupling of the halo dark matter particle. 
For the perturbative case ( A < 0(1)) the above relation implies m ~ lO'^e^. 

To treat particles as a classical field, we require that the inter-particle distance 
should be smaller than their Compton wave length. This gives m ~ 10~'^eV. 

Note that A = Xm^/ AT^m? is very large even for very small A due to the 
smallness of m relative to irip , hence the self-interaction effect is non-negligible. 

Were there any realistic particle physics model satisfying the above relation? 
Unfortunately, the usual cosine potential ^(0) = — Cos(0//)) for the PNGB 
is inappropriate for our study, because the sign of the quartic coupling constant 
is negative in the Taylor expansion about the potential minima and is real. 



Real scalar field such as axion may form oscillating soliton star [III | rather than 
the boson star. 

Instead, we consider the following potential. 

v{<p) = ^,\l + {^ff. (10) 

Inserting mass and quartic coupling from the above potential to the relation in 
eq.d), we get 0.1(^)2 ~ 10^° and equivalently /i ~ lO^el/. 

We also solve the equations numerically and find the dimensionless rotation 



velocity which is given by Vrot = \J M{x)/x = [|(1 — A ^)]2 



The results are shown in fig.l and fig.2. 

Fig.l shows rotation velocity curves for the cases A = and A = 300 . 
The parameters are B{0) = 0.641, (t(0) = 0.1 and B{0) = 0.781, (t(0) = 0.01 
respectively. 

Fig.2 shows a and rotation velocity curve of 8 nodes solution(n = 9). 

It is interesting that the line connecting minimum points of the rotation veloc- 
ity is almost straight. For large n and A ^ 1 mass profile is p ~ r"^ ''. Rotation 
curves are slightly increasing regardless of the self-interactions. 

Including the visible matter may change the slope of the curves and explain 
the variety of the observed galaxy rotation curves as shown in ref.0. 

We will now study the Newtonian limit of our model. The strength of the 
gravity of halo GMhaiol Rhaio is comparable to that of the earth. Therefore we 



can use the Newtonian limit Vt = ^ = y1 + (^)^ 1 , which is comparable 
to the Newtonian gravity approximation 2M{x)/x <C 1. Collecting terms to 




0{^^ = ( — )^) one can find that for A = the equations of motion are[jT5 



V^cr = 7(7, 



(11) 



(12) 



where 7 = 1 — ^. 



Integrating eq.(|l^) and inserting the result into eq. ([TT|) we can find 




(13) 
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which is a dimensionless version of eq.(0). Here E is a integral constant. 

Therefore one may treat the Bose hquid model as a Boson star model. 

It is useful to study the scaling properties of eq. ([TT|) and eq.(|l^) for analyzing 
numerical solutions. 

Rescaling the total number of charges N = Q oc J a'^x^dx I times increases 
the mass / times. 

Equations (0) and (|T2|) are invariant under this rescaling when 



X r^x, a /V, 7 /^7, (14) 

which is consistent with the model in ref.|P]. 

So one may say that for the non-interacting case the heavier halos are smaller 
in size. 

It is difficult to find the scaling properties for the case A 7^ 0. 
For the case A ^ 1 further rescaling : a* = crA^/^, = xA~^/^ and M* = 
MA~^/^ with neglecting terms to 0(A^^) yields wave equations: 

al = (n^B - 1) = -7, (15) 

M: = ^xl{3n'/B + l){nyB-l), (16) 

^' ^ :i-A-') = l-{nyB-iy, (i?) 



ABx, xV 2 



which are also shown in ref.[10 . 



Following the same procedure of A = case, we get the Newtonian limits of 
eq.dH) and eq.(0). 
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whose solutions are 



= 2al = -27, (18) 



Sin{\/2x^) 

7 = -70 7= , (19) 

v2x* 



and 



where 70 = |7(0)|. 

The above approximation is invahd when is large and n > 1. 

As expected, the typical length scale is m~^A5. These solution do not show 
simple scaling property, however numerical study indicates that when the central 
density is lesser than the critical value corresponding to Mmax-, heavier halo has 



smaller radius for both A = and A 7^ cases. [|rj 

Note the facts that the above arguments are valid when the node number is 
fixed and both mass and radius of boson star are increasing function of the node 
number. 



Now, let us discuss the stability of higher node solutions. There are studies 112 



13[| indicating that non-zero node solutions with A = are unstable against fission 
and the small radial perturbation. 

Since higher node solutions are unstable, they must be long-lived to explain 
the age of galaxies. 
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One possible decay mechanism for non-zero node solution to zero node solu- 
tion is gravitaional radiationfl^. The power of gravitaional radiation P is about 



G{cPl /dt^)^, where the quadrapole moment I ~ MhaioR\aio- The available time 
scale is T ~ {R\^i^/ GMhaioY^"^ i which is given by the virial theorem. The param- 
eter a = GMhaiol Rhaio iudicatcs how much halo is relativistic. 

So we can find a crude estimate of power P ~ G^{Mhaio/ RhaioY ~ a^/G ~ 
a^lO^^ergs/ s. Since a ~ 10^'^ for halo, P ~ lO'^^erg/s. The potential energy of 
halo GMl^i^/ Rhaio is about lO^^erg, and therefore time scale of the decay by the 
gravitational radiation is much longer than the age of galaxies ~ lO^^yr ~ lO^^s. 

For the A 7^ case, there is a work indicating that higher node solutions are 



stable against the perturbation with fixed particle number |]19l. However, there 
seems to be no work on the stability against the more general perturbations. 
Since it is still unclear that higher node solution with A 7^ is stable, we must 
again estimate the life time of halo against gravitational radiation. From eq.(15) 
we get Aa^ = —7(7^ , which indicates that in Newtonian limit the energy of the 
repulsive force is comparable to the gravitaional potential energy. So the energy 
distribution in halo is not so different from the A = case. Therefore, we argue 
that the same procedure for calculating the gravitaional radiation is applicable 
to the A 7^= case and halo is long-lived against the gravitaional radiation. 

Another cooling mechanism, evaporation and collapse procedure, is also inefficient! 

In conclusion, we find that self-interactions between the particles, even weak. 
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may play important role in the boson star model of halos. 

Our work can be easily extended to the Boson- Fermion star[jl6| and Q-star 
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Figure Captions 

Figure 1 

Rotation velocities as a function of rescaled x for the parameters A = 0(thin line) 
and A = 300(thick line). Q is 0.9. The real values of the x end are 80 and 220 
respectively. 

Figure 2 

Rotation velocity and 10 x cr as a function of position x ior n — 9 solution. The 

parameters are A = 300, Q = 0.9, B{0) = 0.780 and a{0) = 0.01. 



13 



Figure 1: 
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